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A method is presented by which the attitude of a low-Earth orbiting space- 
craft may be determined using a vector magnetometer, a digital Sun sensor, 
and a mathematical model of the Earth's magnetic field. The method is 
currently being implemented for the Solar Maximum Mission spacecraft (as 
a backup for the failing star trackers) as a way to determine roll gyro 
dri ft. 
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I. INTRODUCTION 


For centuries sailors have used the Earth's magnetic field to guide 
their ships though the oceans of the world by means of the magnetic compass, 
loday it is possible for spacecraft to navigate themselves in much the same 
way, with the mariner's compass replaced by the modern magnetometer. In 
this paper I describe how a vector magnetometer, in conjunction with a 
digital Sun sensor, can be used to determine the attitude of a low-Earth 
orbiting spacecraft. 


This work was motivated by the failure of one of the two star trackers 
on the Solar Maximum Mission (SMM) spacecraft in 1987. The complement of 
attitude sensors on SMM provides only gyroscopes, two star trackers, and 
vector magnetometers for determining the roll attitude. As currently 
written, the on-board computer flight software uses the gyros to determine 
the roll attitude of the spacecraft, with the star tracker used only to 
calculate the roll gyro drift. The magnetometers are not currently used 
for attitude determination on SMM. 

Should the remaining star tracker fail, however, this would leave the 
magnetometers as the only means of determining an absolute roll attitude, 
since the gyros only measure changes in the attitude with respect to 
inertial space. The work described in this paper is a result of an effort 
to determine how SMM's magnetometers might be used as a replacement for the 
remaining star tracker in the event that it fails. 

The approach here will be to find the components of two vectors (the 
cjeomagnetic induction and Sun vectors) in each of two coordinate frames 
(the spacecraft frame and a reference frame); we then solve for the 
rotation matrix between the two frames to determine the spacecraft 
attitude. These calculations will be performed by a computer on the ground 
using data telemetered from the spacecraft; the ground computer will 
calculate roll gyro drift coefficients which will be periodically uplinked 
to the on-board computer. Calculating the gyro drift coefficients on the 
ground will permit ground personnel to select data which was sampled 
while the geomagnetic field was relatively quiet, thus giving the most 
accurate results. 

Section II of this paper describes how the Earth's magnetic field vector 
at the spacecraft position can be calculated from a mathematical model. 
Section III describes the calculation of the Sun vector, and Sections 
IV and V describe how these two vectors may be combined with sensor data 
to determine the spacecraft attitude. Section VI is a summary of the 
paper, and Section VII is a short discussion of associated Leqendre 
functions. 



II. MODELING THE GEOMAGNETIC FIELD 


In order to determine the spacecraft attitude from the magnetometers, 
one must first generate an accurate mathematical model of the Earth's 
magnetic field. Ampere's law at the spacecraft position r is (SI units): 

V « B(r) = 3<r) + (1) 

where ft(r) is the geomagnetic intensity at r, 3(r) is the electric current 
density, and 3^(r)/3t is the displacement current. Since there is no 
current density at r and the geomagnetic field is approximately static. 


we may take 3 and 315/ 3 1 to both be zero. Eq. (1) then becomes 

Vx ft(r) = 3 ( 2 ) 

The constitutive relation for the magnetic induction t is 

l t(r) = Mr) + fi(r) ( 3 ) 

"7 “2 . 

where is the permeability of free space ( 4ir x 10 N A ), and M is the 
magnetization, which is zero at r. Eq. (3) then becomes 

l - Mr) = Mr) ^ 4) 

p. 

Substituting Eq. (4) for ft into Eq. (2) we get 

V x l ft(r) = 15 ( 5 ) 

p. 

since the curl of any gradient is zero, this means that ft can be written as 
the gradient of a magnetic scalar potential V: 1 

Mr) = -y.VV(r) ( 6 ) 


It is conventional in geomagnetism to model the geomagnetic field by 
expanding the magnetic scalar potential V(r) into a Laplace series of 
spherical harmonics with real eigenfunctions: 2 
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V(r,e,A) 


, , . n+l __ _ 

IZ (?) ED 


g nm cos nu + h nm sin mxj P nm (cos 0 ) (7) 


n=l 


m=0 


where 

r is the distance of the spacecraft from the center of the Earth; 

0 is the co-elevation of the sub-satellite point (90* minus the north 
latitude; 

A is the east longitude of the sub-satellite point; 

a is the radius of the Earth, taken to be 6371.2 km; 

P nm (cos e) are the Gauss-normalized associated Legendre functions of 
the first kind; 

g nm and h nm are the Gauss-normalized coefficients of the expansion. 

The n=0 terms in this expansion are absent because they would represent a 
magnetic monopole component of the field; the n=l terms represent the 
dipole component, the n=2 terms represent the quadrupole component, etc. 

The expansion coefficients g nm and h nm are found empirically; they are 

updated every five years and published along with their time derivatives 

(the secular variation) by the International Association of Geomagnetism 

and Aeronomy ( IAGA) . These published coefficients are Schmidt normalized 

and may be used to calculate the geomagnetic scalar potential V(r) using 

Eq. (7) if the Schmidt-normalized associated Legendre functions P m are 

n 

substituted for the Gauss-normalized functions P . Using Gauss normalization 
will save about 7% in computer time, however, ^ so for convenience in 
computer work. Gauss normalization will be used throughout this paper. 

Table I shows the Gauss-normalized coefficients g nm and h nm for the 
International Geomagnetic Reference Field (IGRF) 1985. These were 
calculated from the Schmidt-normalized coefficients published by the 
IAGA and can safely be extended to 1990 with the secular variation 
coefficients in the last two columns. 
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TABLE I. International Geomagnetic Reference Field, IGRF 1985. 
Coefficients are Gauss -no rmal i zed. 




nm 

.nm 

*nm 

*nm 

n 

m 

f - T ) 

( nT ) 

( nT / yr ) 

( nT / yr ) 

1 

0 

-29877 


23.2 


1 

1 

-1903 

5497 

10.0 

- 24.5 

2 

0 

-3110 


- 20.6 


2 

1 

5274 

-3795 

5.9 

- 19.9 

2 

2 

1464 

-268 

6.1 

- 17.5 

3 

0 

3250 


12.8 


3 

1 

-6761 

-955 

- 14.1 

16.2 

3 

2 

2409 

550 

- 1.2 

4.5 

3 

3 

660 

-234 

0.1 

- 8.5 

4 

0 

4099 


0.4 


4 

1 

4317 

1289 

- 3.3 

21.0 

4 

2 

1420 

-978 

- 30.5 

8.6 

4 

3 

-891 

142 

- 2.9 

5.2 

4 

4 

125 

-220 

- 5.0 

0.7 

5 

0 

-1693 


10.2 


5 

1 

3619 

478 

1.0 

1.0 

5 

2 

1944 

1137 

- 11.5 

- 1.5 

5 

3 

-442 

-729 

- 15.1 

- 0.5 

5 

4 

-357 

-166 

0.2 

1.3 

5 

5 

-34 

67 

- 0.1 

0.0 

6 

0 

751 


20.2 


6 

1 

1229 

-302 

- 5.7 

- 7.6 

6 

2 

747 

1345 

25.4 

- 16.4 

6 

3 

-1853 

687 

6.0 

- 8.0 

6 

4 

22 

-273 

0.0 

- 12.6 

6 

5 

40 

-9 

2.1 

- 1.2 

6 

6 

-69 

13 

0.8 

- 0.1 

7 

0 

2011 


5.4 


7 

1 

-2164 

-2909 

- 21.3 

7.1 

7 

2 

58 

-753 

- 14.5 

29.0 

7 

3 

491 

-20 

16.4 

22.5 

7 

4 

-74 

284 

12.3 

23.5 

7 

5 

25 

105 

2.5 

1.9 

7 

6 

22 

-51 

- 1.2 

0.5 

7 

7 

0 

-4 

- 0.1 

0.6 

8 

0 

1056 


35.2 


8 

1 

402 

469 

0.0 

6.7 

8 

2 

0 

-1178 

16.8 

- 56.1 

8 

3 

-456 

207 

16.6 

4.1 

8 

4 

-241 

-668 

- 8.0 

- 21.4 

8 

5 

30 

163 

- 4.4 

3.0 

8 

6 

27 

82 

0.7 

- 5.5 

8 

7 

10 

-40 

- 1.3 

- 0.3 

8 

8 

-4 

-6 

- 0.5 

0.8 
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TABLE I (cont.) 




nm 

q 

h nm 

n 

m 

(nT) 

(nT) 

9 

0 

475 


9 

1 

1274 

-2675 

9 

2 

109 

1738 

9 

3 

-996 

747 

9 

4 

507 

-282 

9 

5 

-101 

-202 

9 

6 

-17 

157 

9 

7 

53 

75 

9 

8 

5 

-16 

9 

9 

-3 

1 

10 

0 

-722 


10 

1 

-973 

243 

10 

2 

421 

0 

10 

3 

-826 

496 

10 

4 

-234 

701 

10 

5 

370 

-296 

10 

6 

124 

0 

10 

7 

20 

-20 

10 

8 

16 

33 

10 

9 

8 

0 

10 

10 

0 

-4 


(nT/yr) 


h nm 

(nT/yr) 
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When the coefficients in Table I are substituted into Eq. (7), we get 
the geomagnetic scalar potential V(r); substituting this V(r) into Eq. (6) 
yields the geomagnetic induction vector ^(r) in a reference frame fixed 
in the Earth (which will be referred to as the EB (Earth-based) frame). 

The EB frame has its origin at the center of the Earth, its x axis pointing 
out of the intersection of the equator with the prime meridian, its z axis 
pointing out of the Earth's geographic north pole, and its y axis in the 

/\ A 

x x z direction. 


If we calculate the gradient in Eq. (6) in spherical polar coordinates, 

m = a e + 1 21 e + —4— 21 S, (8) 

Vv e r r 99 0 r sin 0 9X A 


the resulting spherical components of B may be used to easily calculate the 
standard geomagnetic elements: 2 


X = -Bq 

Y = b a 

Z = -B 

H = (X 2 + Y 2 )** 

V = |z[ 

F = (X 2 + Y 2 + l 2 ) h 

D = arctan (Y/X) 

I = arctan (Z/H) 


(Northward component) 
(Eastward component) 
(Downward component) 
(Horizontal intensity) 
(Vertical intensity) 
(Total intensity) 
(Magnetic declination) 
(Magnetic inclination) 


A computer program which calculates 1i(r) can then be checked by comparing 
the geomagnetic elements it calculates with the elements found in charts 
and tables in the literature. 5 


For spacecraft attitude determination, we will need to know the 
components of the modeled geomagnetic induction vector l(r) in the 
geocentric inertial (GCI) reference frame rather than the EB frame. 

The GCI frame is fixed with respect to the stars and has its origin at the 
Earth's center, its x axis toward the vernal equinox, its z axis out of 
the Earth's geographic north pole, and its y axis in the z x x direction. 
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The GCI frame differs from the EB frame only by a rotation about their 
conmon z axis. Specifically, in cartesian coordinates 

^GCI = R B EB (9) 

where the rotation matrix R is given by 


R 


cos y sin y 0 

-sin y cos y 0 

0 0 1 


( 10 ) 


and where y is the Greenwich hour angle of the vernal equinox, which is 
equal to the sidereal time at Greenwich (GST) and is given by 6 

Y - LST - (X/15* ) (11) 

where LST is the local sidereal time and X the east longitude of any 
convenient point on the Earth's surface. An expression for y which is 
often more convenient Is 6 

Y = 99*. 6910 + 36000*. 7689T + 0* .0004T 2 + UTC (12) 

where T is the time (in Julian centuries of 36525 days) since 1900 and 
UTC is Coordinated Universal Time expressed in degrees. 

Alternatively, if we work in spherical polar coordinates, Eqs. (9) and 
(10) may be replaced with 


GCI 

r EB 


GCI 

= 

6 EB 

(13) 

GCI 

= 5 - Y 

a eb 
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III. SUN VECTOR CALCULATION 


The geomagnetic induction vector ff(r) modeled in Section II is by itself 
insufficient to determine the spacecraft attitude, since if the spacecraft 
is rotated about the t vector it will still yield the same components of t 
in a reference frame fixed in the spacecraft; hence one degree of freedom 
is left unspecified. 7 It is therefore necessary to know the components of 
one more vector (not parallel to §(r)) in order to specify the spacecraft 
attitude completely. For this discussion we choose the Earth-to-Sun vector 
(or simply the "Sun vector") S, and in this section I discuss how to 
calculate the components of S in the GCI frame. (The cartesian GCI compon- 
ents of S will be used along with the cartesian GCI components of B from 
Section II to determine the spacecraft attitude in Sections IV and V.) The 
Sun unit vector S is given in cartesian GCI coordinates approximately by 
(ignoring the small corrections for parallax and light aberration)® 


S(t) 


COS L 

sin L cos e 
sin L sin e 


(14) 


where L{t) is the mean longitude of the Sun and e is the mean inclination 
of the ecliptic from the Earth's equatorial plane. 

The mean longitude of the Sun l{t) may be calculated from 8 

L(t) = L{t ) + M{t) + 2e sin M{t) + 6 1 U 5 ) 


where t is a reference time, e is the eccentricity of the Earth's orbit 
(e = 0.016722), M{t) is the mean anomaly of the Sun, and B is defined by 


B 




(16) 


where t is the length of the sidereal year (3.1558149548 x 10 seconds) 

and t is the length of the anomalistic year (3.1558433 x 10 seconds), 
d y 

The mean anomaly of the Sun M[t) is given by 8 


M{t) 


M[t ) 
0 



a) 


(17) 
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where o is the time it takes light to travel from the Sun to the Earth, 
about 499 seconds. 

Finally, the mean inclination of the ecliptic e is given by 8 
e = 23*27' 8.26" - 46.845" T ( 

where T is the time in Julian centuries since 1900. 


IV. GENERAL ATTITUDE DETERMINATION 

^ Knowing the cartesian components of the geomagnetic induction unit vector 
B and the Sun unit vector S in the GCI reference frame (as described in 
Sections II and III), we can now determine the spacecraft attitude if we 
also know the components of these two vectors in the spacecraft (SC) 
reference frame (i.e. along the roll (a?) , pitch (y) , and yaw (z) axes). 

The essence of the problem of determining the spacecraft attitude is then 
this: we are given the components of two vectors in each of these two 

reference frames which have a common origin, and we must solve for the 
rotation matrix between the two frames. There are several ways of doing 
this; I will describe here one of the simplest methods, known as the 
algebraic method . 8 Let £ be some vector whose components we wish to 
transform from the GCI reference frame to the SC frame. The two sets of 
components are related by 

*£■ _ - 

4 SC " A 5 GCI (19) 

where | $c and t GCI are column vectors containing the SC and GCI components, 
respectively, and A is the rotation matrix we wish to solve for. Formally , 
we could solve this equation of A by post-multiplying both sides by the 
inverse of the matrix 

G LI 

-K + “I 

A " 5 SC ? GCI 
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Unfortunately, c QCI is not a square matrix, so we cannot take its inverse 
directly (without a somewhat messy diversion into pseudo-inverses). 



In the algebraic method, we use the two known vectors B and S to 
construct an orthogonal triad of vectors (a, b, and c) in each reference 
frame (SC and GCI). (This will work only if B and S are not parallel.) 

We then define a matrix M for each frame which has the vectors a, b, and c 
as its columns; the columns of this matrix will transform from GCI to SC 
coordinates just like the vector t above, and we will be able to take its 
inverse (since it will be a square matrix), thus allowing us to solve for 
the rotation matrix A. 

Specifically, let us define the vectors a, b, and c in the GCI and SC 
reference frames as follows: 


GCI Frame SC Frame 


GCI 

S GCI 

a SC 

b SC 


^GCI X ®GCI 

^SC 

^SC x ^sc 

'GCI 

I^GCI X ®GCI 1 

l^sc x ^SC 1 

GCI 

a GCI X ^GCI 

/\ 

c SC 

*sc x *sc 


where : 

B is the geomagnetic induction unit vector in GCI cartesian 
GCI 

coordinates, calculated as in Section II# 

S is the Sun unit vector in GCI cartesian coordinates, calculated 
GCI 

as in Section III. 

B is the geomagnetic induction unit vector in SC cartesian 
coordinates, from the spacecraft magnetometers. 

S is the Sun unit vector in SC cartesian coordinates, from the 
S C 

spacecraft digital Sun sensor. 
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We now construct two 3x3 square matrices M QCI and M $c whose columns are 
the cartesian components of a, b, and c: 

M„‘ t = r a... I b„_ I r 1 


Since the rotation matrix A rotates column vectors from the GCI to the SC 

reference frame, it will also rotate each column of M rrT into the corre- 
. btl 

sponding column of M sc : 


M SC ‘ A M GCI 


Since is a square matrix, we can now solve for the rotation matrix A 

by post-multiplying both sides by the inverse of M„~ x : 

GCI * 

A = M sc m gci 1 


Furthermore, since Mggj was defined to be an orthogonal matrix, its inverse 
is equal to its transpose and so 


A 



( 22 ) 


The matrix A given by Eq. (22) rotates any vector from the GCI to the 
spacecraft reference frame and thus determines the spacecraft attitude. 
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V. SMM ROLL ATTITUDE DETERMINATION 


On SMM, the roll, pitch, and yaw angles are measured with respect to 
the somewhat whimsically named "SUN" reference frame, defined by the three 

A A A 

orthonormal vectors S, U, and N: 


S is the Sun unit vector, calculated as in Section III. 

* _ p x S (0 points out of the east limb of the Sun, (23) 

u = 7* *7 along the solar equator.) 

1 P X s ] 

N=SxU (NisP projected onto the plane normal to the 

line of sight.) 

where P is the solar spin axis unit vector. In order to determine the 
spacecraft attitude from the magnetometers and the Fine Pointing Sun 
Sensor (FPSS), we will need to know the components of the vectors S, U, 
and N in the GCI reference frame. 

We first need to calculate the cartesian components of the solar spin 
axis unit vector P in the GCI frame. We begin by working in the eoliptio 
(ECL) frame, defined by the unit vectors * ECL . J/ E CL » and 3 ECL : 

av r , points in the direction of the vernal equinox; 

LLL a a 

y ECL points in the direction of 3 £CL x ® ECL ; 

2 ECL points toward the ecliptic north pole. 

The vector P in the ecliptic frame has cartesian components 8 


P 


ECL 


s i n s i n i 
-cos Q. sin i 
cos i 


(24) 


where ft is the longitude of the ascending node of the solar equator on the 
ecliptic, given by 8 


ft = 73* 40' + 50* 25' A 


(25) 
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where A is the time in years since 1850; i is the inclination of the solar 
equator to the ecliptic: 


i = T 15' 


(26) 


The ecliptic and GCI reference frames differ only by a rotation of 
magnitude -e about their common x axis, where e is the inclination of the 
ecliptic (Eq. (18)). Hence the GCI cartesian components of P are 


P 

GCI 


1 0 0 

0 cos e -sin e 

0 sin e cos e 

J 



(27) 


Having found the cartesian components of the solar spin axis unit vector P 
in the GCI frame, and knowing the components of S from section III, we may 
now use Eq. (23) to determine the cartesian components of the U and N 
vectors in the GCI frame; these will be used in Eq. (29) below to determine 
the SMM roll attitude. 

It is especially convenient to determine SMM's roll angle p when the 
spacecraft's FPSS is pointed at Sun center (so that the pitch and yaw 
angles are zero). In this case the x axis of the spacecraft reference 
frame (the roll axis) will coincide with the $ vector of the SUN frame. 

The two reference frames will then differ only by a rotation about their 
common x-S axis (Fig. 1). The roll angle p is then given by 8 


P = 4> i - 4> 2 (28) 

where cp ^ is calculated from the magnetic field model and solar ephemeris 
(Sections 1 1 -V) : 




l 


•i 

tan 


B GCI * U GCI 
B GCI * ^GCI 


(29) 
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and <t> ^ is found from magnetometer data: 


<P 

2 


■1 


tan 


mag 


pitch 


1 


mag 


yaw 


(30) 


Substituting Eqs. (29) and (30) into Eq. (28) then determines the SMM roll 
attitude. Note that we do not require any FPSS data in this case; the Sun 
vector S^ was tacitly assumed to lie along the S axis of the SUN frame 
since the FPSS is pointed at Sun center. 


VI. SUMMARY 

Two methods for spacecraft attitude determination using the Earth's 
magnetic field have been presented: the algebraic method for spacecraft 

in general (Section IV), and a simpler method specific to SMM (Section V). 
Both methods compare the Earth's magnetic field as calculated by a 
mathematical model (Eq. (7)) with magnetometer measurements; in addition, 
the Sun vector calculated by Eq. (14) and measured with a digital Sun 
sensor is used to completely specify the spacecraft attitude. 


VII. APPENDIX: ASSOCIATED LEGENDRE FUNCTIONS 


The associated Legendre functions of the first kind form a complete 
orthogonal set of functions over the interval 0 = [0, tt! ; it is this 
property which makes them a useful basis in which to expand the geomagnetic 
scalar potential V(r). 

Although a variety of normalization conventions for the associated 
Legendre functions are in common use, the three most common are the so- 
called Neumann, Schmidt, and Gauss normalization conventions. Neumann 
normalization is the convention most often found in mathematics text- 
books! 0-12 The geomagnetic coefficients g and h given in the literature 
are usually defined for the Schmidt-normalized associated Legendre 
functions; Schmidt normalization has the advantage that the normalization 
constants are independent of m (for m f 0) for any given n, so the 
relative strengths of the different terms can be easily judged. Gauss 
normalization is useful because it saves about 1% in computation time 
on a computer 3 ; the g and h coefficients in this paper have been converted 
to Gauss normalization for ready use in computer work. 

The various normalization conventions are defined as follows: 


Neumann normalization (P ) 


/. 


P (cos 0 ) P„ (cos e) sin e de 
nm £m 


Schmidt normalization (P p m ) 


/. 


P m (cos e) P m (cos e) sin 0 d© 

n ic 


Gauss normalization (P nm ) 


f 


P nm (cos 0) P x, "(cos 0) sin e d0 


Jm 


2 (n+m)! * 

2n + 1 (n-m) ! n£ 


2 ( 2 - 5 m0> 

"TrT+1 — 


6 


n£ 


2 (n-m) ! (n+m) ! r 

(2n+l)I('2n-TTrTF n£ 
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One may easily calculate the conversion factors for converting between 

Schmidt and Gauss normalizations by simply taking the square root of the 

quotient of the respective normalization constants. In particular, if 

the conversion factors S are defined by 13 

nm J 


then 


2 ( 2 - 


‘W 


nm 


*5 


~'~2n + 1 °n£ 

2(n-m)! (n+m)! " 

(7n+l) [(2n-l) ! TP 6 n* 


or 


S 


nm 


~(2-6m 0 ) ("-">)! 

(n+m) ! 



(2n-l ) I ! 
"Ih-mJ! 


Table II lists the explicit values of these conversion factors up to n=12. 

To convert a table of Schmidt-normalized coefficients g m and h m (such as 
. 3 n n v 

those usually found in the literature) from Schmidt to Gauss normalization 

use 


nm 

9 

h nm 



m 

m 


and analogous expressions for the secular variation coefficients. 



TABLE II. Conversion factors S between Schmidt- and Gauss- normalized 

nm 

associated Legendre functions of the first kind. 


n 

m 

S 

n 

m 

S 



nm 

— 

— 

nm 

1 

0 

1 

9 

0 

12155/128 

1 

1 

1 

9 

1 

(7293/128) /5 

2 

0 

3/2 

9 

2 

(663/64) /110 

2 

1 

/3 

9 

3 

(221/128) /2310 

2 

2 

(1/2) /3 

9 

4 

(51/64) /5005 

3 

0 

5/2 

9 

5 

(255/128) /286 

3 

1 

(5/4) /6 

9 

6 

(17/64) /4290 

3 

2 

(1/2) AS 

9 

7 

(51/256) /1430 

3 

3 

(1/4) /10 

9 

8 

(3/128) /1 21 5 5 

4 

0 

35/8 

9 

9 

(1/256) /24310 

4 

1 

(7/4) /10 

10 

0 

46189/256 

4 

2 

(7/4) /5 

10 

1 

(4199/128) /55 

4 

3 

(1/4) /70 

10 

2 

(4199/256) /165 

4 

4 

(1/8) /35 

10 

3 

(323/128) /4290 

5 

0 

63/8 

10 

4 

(323/128) /2145 

5 

1 

(21/8) /15 

10 

5 

(323/128) /858 

5 

2 

(3/4) /105 

10 

6 

(323/512) /4290 

5 

3 

(9/16) /70 

10 

7 

(19/256) /72930 

5 

4 

(3/8) /35 

10 

8 

(19/256) /I 2155 

5 

5 

(3/16) /14 

10 

9 

(1/256) /461890 

6 

0 

231/16 

10 

10 

(1/512) Z92378 

6 

1 

(33/8) /21 

11 

0 

88179/256 

6 

2 

(33/32) /210 

11 

1 

(29393/512) /66 

6 

3 

(11/16) /210 

11 

2 

(2261/256) /2145 

6 

4 

(33/16) /7 

11 

3 

(969/512) /30030 

6 

5 

(3/16) /154 

11 

4 

(969/128) /1 001 

6 

6 

(1/32) /462 

11 

5 

(6783/512) /143 

7 

0 

429/16 

11 

6 

(399/512) /14586 

7 

1 

(429/32) /7 

11 

7 

(133/512) Z36465 

7 

2 

(143/32) /42 

11 

8 

(7/256) Z692835 

7 

3 

(143/32) /21 

11 

9 

(21/512) /46189 

7 

4 

(13/16) /231 

11 

10 

(1/512) /I 939938 

7 

5 

(13/32) /231 

11 

11 

(1/512) /88179 

7 

6 

(1/32) /6006 

12 

0 

676039/1024 

7 

7 

(1/32) /429 

12 

1 

(52003/512) /78 

8 

0 

6435/128 

12 

2 

(7429/512) /3003 

8 

1 

2145/32 

12 

3 

(7429/512) /2002 

8 

2 

(429/64) /70 

12 

4 

(22287/2048) /2002 

8 

3 

(39/32) /II 55 

12 

5 

(1311/512) /17017 

8 

4 

(195/64) /77 

12 

6 

(3059/1024) /4862 

8 

5 

(15/32) /1001 

12 

7 

(161/512) /1 38567 

8 

6 

(15/64) /858 

12 

8 

(161/1024) /I 38567 

8 

7 

(3/32) /715 

12 

9 

(23/512) Z323323 

8 

8 

(3/128) /715 

12 

10 

(23/1024) /1 7 6358 
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11 

(1/512) /2028117 




12 

12 

(1/2048) /1352078 
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